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Abstract. The distribution of TeV spectral slopes versus redshift for currently known TeV blazars (16 sources with z<0.21, 
and one with z>0.25) is essentially a scatter plot with hardly any hint of a global trend. We suggest that this is the outcome 
of two combined effects of intergalactic 77 absorption, plus an inherent feature of the SSC (synchro-self-Compton) process 
of blazar emission. First, flux dimming introduces a bias that favors detection of progressively more flaring sources at higher 
redshifts. According to mainstream SSC models, more flaring source states imply sources with flatter TeV slopes. This results in 
a structured relation between intrinsic TeV slope and redshift. The second effect, spectral steepening by intergalactic absorption, 
affects sources progressively with distance and effectively wipes out the intrinsic slope-redshift correlation. 
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, 1. Introduction 



With the increasing number of very high energy (VHE) 7- 
ray detection of blazars, a broad correlation between observed 
slope and redshift was generally expected, in the sense that 
higher redshift sources would tend to show systematically 
steeper spectra. 

The reason for this is absorption of the emitted VHE emis- 
sion by the intergalactic background light (IBL), i.e. the in- 
tegrated optical/IR emission from the evolving stellar popu- 
lations in galaxies (e.g., Hauser & Dwek 2001). In fact, pair 
creation due to the interaction of a hard photon (with energy 
E) with a soft photon (with energy e) is expected provided 



The cross section of the 77— >e 



interaction, 

-2\2 , 



that eE>m e c 

a 7J (E, e) (Heitler 1960) is maximized when e~2(m e c 2 ) 2 / E. 
Then, if £~1 TeV, pair creation will be most likely for 
e^0.5 (E /TeV) -1 eV, which corresponds to a 2 /zm (A"-band) 
photon. The formalism required to evaluate this effect has been 
largely developed by Stecker and collaborators (e.g., Stecker 
1971; Stecker et al. 1992 and 2006). 

However, when we plot the currently known observed TeV 
slopes against redshift (Fig. 1), we see a scatter plot with hardly 
any sign of a global correlation: neglecting slope uncertainties, 
a global linear correlation is suggested only marginally (~2cr) 
by the whole data sample (r=0.52, N=18). Even so, most of the 
suggestion's strength comes from only the two points at z>0: 
for the bulk of the distribution, which provides a fair sampling 
of the redshift interval z<0.2, there is no correlation (r=0.22, 
N=16). In this note we propose an interpretation of Fig. 1. 
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Fig. 1. Observed TeV slope vs redshift for known TeV blazars. 
Bars denote statistical uncertainties as quoted in the original 
papers. 



2. Intergalactic absorption 

The cross section for the reaction 77^e ± is (Heitler 1960), 
er 77 (£,e) = 1.25 x 10~ 25 (1 - (3 2 ) x 



x 2/3(/3 2 -2) + (3-/3 4 )ln 



1-/3 



(1) 
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where /3=a/1 — (m e c 2 ) 2 / Ee. By calling n(e) the number per 
unit volume of background photons with energy equal to e at 
the current position, the corresponding optical depth due to the 
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Table 1. TeV blazar data. 
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Col.l: source name. 
Col.2: source redshift. 

Col.3: observed 0.2-2 TeV photon spectral index, and associated statistical uncertainty. The corresponding systematic uncertainties are typically 
~0.1 for H.E.S.S. and ~0.2 for MAGIC. 

Col.4: observed >0.2 TeV flux (in erg cm~ 2 s -1 ), and associated statistical uncertainty (from observed spectral normalization only). 

Col.5: Cherenkov telescope (CT) or array with which the data in col.3 have been collected: symbols stand for H=H.E.S.S., M=MAGIC, 

W=Whipple. 

Col.6: corrected 0.2-2 TeV photon spectral index. 
Col.7: corrected >0.2 TeV flux. 

Col.8: 5 GHz flux density (in Jy): data are from the NED, except for 1ES 0347-121 (Fossati et al. 1998) and H 2356-309 (Costamante & 
Ghisellini 2002). 

Col.9: 2-10 keV flux (in erg cm~ 2 s _1 ): the data come from the papers listed in Colli or from references quoted therein; when necessary, 
fluxes have been converted from A'XTE/ASM count rates using the conversion 1 ct/s =3.33E— 10 erg cm~ 2 s _1 . 
Col.10: minimum 2-10 keV flux (in erg cm~ 2 s _1 ). 

Colli: Reference for VHE 7-ray data, and for the baseline X-ray flux when the latter was not taken from the NED - Al+07a: Albert et al. 
2007a; Al+07b: Albert et al. 2007b; Al+07c: Albert et al. 2007c; Al+06a: Albert et al. 2006a; Al+06b: Albert et al. 2006b; Al+07d: Albert et 
al. 2007d; Ah+05a: Aharonian et al. 2005a; P+99: Perlman et al. 1999; Ah+08: Aharonian et al. 2008; D+01: Donato et al. 2001; Ah+05b: 
Aharonian et al. 2005b; Hor+02: Horan et al. 2002; S+97: Sambruna et al. 1997; Ah+07c: Aharonian et al. 2007c; P+96: Perlman et al. 1996; 
D+05: Donato et al. 2005; Ah+06a: Aharonian et al. 2006a; W+84: Wood et al. 1984; Al+06c: Albert et al. 2006c; Ah+07a: Aharonian et al. 
2007a; Ah+07b: Aharonian et al. 2007b; Al+07e: Albert et al. 2007e; Al+07f: Albert et al. 2007f. 

attenuation between the source redshift, z e , and the Earth, is 

c r z - f 2 x 

t 17 {E) = — / s/1 ± z dz / -dx x 

-"0 Jo JO 1 

/>oc 

x / n(e) o-, n (2xEe(l + z) 2 ) de , (2) 

J 2 (m e c ) 

where x=(l— cos8), 9 being the angle between the pho- 
tons, and 7J is the Hubble constant. Eq.(l) assumes 17 = 1 



and no redshift evolution of n(e) - the latter assumption 
being adequate within the relatively low redshifts relevant 
to this paper. For demonstration purposes let us assume, 
following Stecker et al. (1992), that the local IBL has a 
power-law form, n(e)oce -2 ' 55 : then the above integral yields 
t(E, z)cx^ 1 - 55 z^ with ?7-1.5 (see Stecker et al. 1992). This 
calculation shows an important property of r 77 : it depends 
both on the distance traveled by the hard photon (hence on 
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z) and on its energy.The expected VHE 7-ray flux at Earth 
will be, of course: F(E)= (d//dE) e^rr^) (differential) and 
F{>E)= J™(dI/dE') e- r -rr( E 'UE' (integral). 

Along the same path, but using recent data on galaxy lu- 
minosity functions and redshift evolution, Stecker et al. (2006, 
2007) made a detailed evaluation of the IBL density as a func- 
tion of both energy and redshift for 0.003 eV < e <13.6 eV 
and 0<z<6 in the ACDM universe with f2A=0.7 and il, n =0.3. 
Using their calculated IBL photon densities, they calculated 
r 77 for 7-rays with energies 0.004<(£YTeV)<100 emitted by 
sources at redshift 0<z<5. 

Stecker & Scully (2006) fitted such t 11 (E 1 z) to a form 
that was assumed to be logarithmic in E in the energy range 
0.2 TeV<£X2 TeV and linear in z for 0.05<z<0.4, i.e., 
T yy (E, z)={A + Bz) + (C + Dz) \n(E/TeY). Consequently, 
if the intrinsic source spectrum can be described as a 
power law, F s (E)—KE^ a , in the band between 0.2 TeV 
and 2 TeV, the observed emission can be approximated as 
F {E)=Ke-( A+B ^E- < - a+c+Dz \ i.e., the spectrum will be 
attenuated by a factor e -( A + Bz ) and steepened by Aa—C + 
Dz. The numerical values given by Stecker & Scully (2006) are 
A= -0.346 (-0.475), B=16.3 (21.6), C= -0.0675 (-0.0972), 
and D=7.99 (10.6) in the Stecker et al. (2006) baseline (fast) 
evolution model. 

In Fig.2a,b we plot the IBL-corrected slopes versus red- 
shift, using the baseline (Fig. 2a) and fast (Fig. 2b) evolution 
models of Stecker & Scully (2006) Q. A loose trend of flatter 
slopes at higher redshift is apparent now: from Table 1 we de- 
rive <a 7 orr >=2.62±0.16forz<0.1 and <a 7 orr >=1.97±0.15 
for z>0.1 (using the baseline correction; the fast-evolution cor- 
rection only enhances the effect). In this paper we suggest that 
the trend of intrinsic slope with redshift is induced by a lumi- 
nosity effect - the latter, in turn, resulting from a selection ef- 
fect. A possible correlation of spectral flatness and luminosity 
with redshift was also pointed out by Stecker et al. (2007). 

3. A luminosity effect 

Let us consider which emission conditions in blazars can lead 
to flatter TeV slopes. The standard synchrotron-Compton emis- 
sion model of blazars (e.g., Jones et al. 1974; Maraschi et al. 
1992; Mastikiadis & Kirk 1997) gives us a clue on the link be- 
tween active phase and TeV spectral shape in blazars. Its basic 
feature is that, if the distribution function of the emitting elec- 
trons changes at the highest energies, correlated flaring at X-ray 
and TeV energies are produced, with the highest energy elec- 
trons producing X-rays via snchrotron and the TeV radiation 
via IC scattering. 

A well-known example of this correlation is provided by 
Mrk 501 (Pian et al. 1998): during the giant flare of April 
1997, its (apparent) bolometric luminosity increased by a fac- 
tor of ^20, its synchrotron and IC peaks shifted to higher en- 
ergies by, respectively, a factor ^100 (into the hard X-ray 

1 Mrk 421, Mrk 501, and PKS 2344+514 are all at z<0.05, which 
is the lower limit for the analytic approximation given in Stecker 
& Scully (2006): this, although it means slightly overestimating the 
(modest) true absorption for these sources, will in no way affect the 
main conclusions of this paper. 



range, £ 50-100 keV) and a factor 10 (Klein-Nishina lim- 
ited, into the sub-TeV range), and /2-iokcV increased by a 
factor of ^3. As the shift of the synchrotron peak by ^2 or- 
ders of magnitude could not be due to variations of the mag- 
netic field and/or the Doppler factor because enormous varia- 
tions would have been demanded on either quantity that would 
have revealed themselves by affecting other parts of the spec- 
trum (which was not observed), the luminosity and spectral 
variation could be explained by invoking an injection, continu- 
ous throughout the duration of the burst, of energetic electrons 
that became superposed to the baseline distribution of emitting 
electrons responsible for the quiescent emission of Mrk 501. 
The former had a spectrum that was flatter than the latter, and 
extended to higher (by a factor ~ 10-30) energies (see Pian et 
al. 1998). The resulting spectral energy distribution (SED) in- 
volved, w.r.t. the baseline case, a higher power output, a higher 
Compton/synchrotron power ratio, and both the synchrotron 
and Compton peaks located at higher energies: that meant a 
harder TeV slope and enhanced 2-10 keV and TeV fluxes. The 
hardening and increase of the TeV emission, simultaneously 
with an increase of the 2-10 keV emission, has been observed 
also in other blazars, down to very short (e.g., hours and min- 
utes) time scales (e.g., see Mrk421 and Mrk501 in Table 1 and 
references cited therein). 

Because all the known TeV blazars, including Mrk 501, 
are of the HBL type (with one LBL exception, BL Lacertae), 
whose SEDs are well reproduced by the SSC model (Ghisellini 
et al. 1998), we expect that the basic features of Pian et al.'s 
(1998) model of the 1997 burst of Mrk 501 apply to TeV blazar 
in general. Hence, we do expect a correlation between harder 
TeV spectral shapes and more active emission states in (HBL) 
blazars. We then have to examine whether more distant sources 
are found, on the average, in more active states. 

4. A distance bias for active states? 

Within the theoretical SSC framework outlined above, two 
indicators of enhanced emission activity are (i) the TeV/X- 
ray luminosity ratio, and (ii) the X-ray enhancement. The 
former quantity, based on simultaneous data, represents the 
Compton/synchrotron power ratio. The latter quantity is de- 
fined as the ratio of the 2-10 keV power, measured simul- 
taneously with the TeV observations, to its lowest historical 
level: it measures the increase of the synchrotron power, i.e. 
the strength of the activity. 

( i) In Fig. 2c we plot the ratio of intrinsic TeV to X-ray lumi- 
nosity versus redshift. In both bands, the luminosities are com- 
puted from the observed fluxes and redshifts (k-corrections do 
not significantly alter the results ), whereas the TeV lumi- 
nosities are IBL-corrected using the baseline model of Stecker 
& Scully (2006). Even though the fields of view (FoVs) of the 

2 The full-band k-correction, to be applied to the 2-10 keV lumi- 
nosities and the 0.2-2 TeV luminosities, amounts to a factor (l + z) r ~ 2 
if in the relavant band the spectrum is power law with intrinsic slope 
F. For the blazars in Table 1, it is either z<0. 1 or a 7 orr ~2 (and typ- 
ically Qf2-iokcV~2.2, see Donato et al. 2005), so the corresponding 
k-corrections to the X-ray and VHE 7-ray luminosities are small. 
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Fig. 2. Parameter trends for the TeV blazars in Table 1 . (i) IBL-corrected VHE slope vs redshift. The spectral slopes have been 
corrected for IBL absorption according to the 'fast evolution' (a) and 'baseline' (b) IBL models of Stecker et al. (2006) as 
parameterized by Stecker & Scully (2006). (ii) VHE to 2-10 keV luminosity ratio versus redshift (c), and 2-10 keV luminosity 
versus redshift (d). Given the modest redshifts encompassed by the current sample of TeV blazars, distances have been computed 
according to D(z)=cz/Hq (with H =72 km s" 1 Mpc" 1 ), and luminosities have not been k-corrected. The VHE luminosities 
are IBL-corrected using the baseline model of Stecker & Scully (2006). In panel (d), the filled dots represent X-ray data taken 
simultaneously with the TeV measurements, whereas the empty squares represent historical low flux states; the thin dotted curve 
represents the 'visibility function' for a 2-10 keV instrumental sensitivity of 7xl0" 12 erg cm" 2 s" 1 (<~0.3 mCrab). (Hi) TeV 
luminosity versus redshift (e) and versus TeV spectral slope (f). Luminosities and slopes are IBL-corrected using the baseline 
model of Stecker & Scully (2006). (iv) Intrinsic (g) and observed (h) TeV slopes versus 5 GHz monochromatic radio luminosities. 



different X-ray and Cherenkov detectors are different, since the 
measured blazar fluxes do originate from (effectively) point- 
like sources and no other (known) X-ray or TeV sources lie 
in the detectors' FoVs, in Fig. 2c there is no need to account 
for the different FoVs. Furthermore, in the SSC framework the 
synchrotron and Compton emissions are essentially cospatial 
and hence are affected by the same beaming-related boost fac- 
tor. Indeed there is a sharp correlation, suggesting an increasing 
importance of the Compton hump for more distant sources. 

( ii) In Fig. 2d we plot the X-ray luminosities of our TeV blazars 
versus redshift (filled dots: fluxes measured simultaneously 
with the TeV measurements; open squares: faintest historical 
fluxes). The lower limit to the allowed L2-iokev(z) distri- 
bution in Fig. 2d is dictated by a combination of X-ray tele- 
scope sensitivities, exposure times, and source spectral shape: 
for sake of illustration, in Fig. 2d we plot an effective limiting 
sensitivity of / s =0.7x 10" 11 erg cm" 2 s" 1 . A trend is apparent 
in the data. Whereas the ground-state emissions (open squares) 
tend to be found close to the limiting curve, the TeV-selected 
X-ray emissions tend to lie significantly above it. The strongest 
enhancements, as defined here, occur for the nearest sources 
(z<0.1): this however stems from a selection effect, as fainter 
fluxes can be detected from more nearby, but otherwise simi- 
lar, sources. For more substantial redshifts (z>0.1), we notice 
that the X-ray emissions tend to be more enhanced w.r.t. their 
ground states with increasing z: although noisy, there appears 
to be a link between such enhancements and redshift. 



Thus, both our activity indicators suggest that TeV sources 
further away tend to be observed during more active phases. 
Why? Our tentative answer - a selection effect - is based on 
the following argument. Moving to higher redshifts, X-ray and 
TeV measurements are both affected by geometrical flux dilu- 
tion (oc z" 2 ), but the latter are also affected by IBL aborption 
(oc e~ T i'<^). This implies that, in order for sources to be de- 
tected at both ke V and TeV energies with a given pair of X-ray 
and Cherenkov detectors, the ratio of intrinsic TeV to 2-10 keV 
luminosities must increase with redshift starting from a thresh- 
old redshift that depends on the sensitivity of the two instru- 
ments. Recalling the behavior of the SSC model for Mrk 501's 
giant 1997 flare by Pian et al. (1998), this means observing 
sources that are in increasingly higher states of emission. 

This argument entails a testable prediction. If distant TeV 
blazars are observed during strong flares, then, because of the 
high electron energies involved, the electrons' scattering cross- 
section is reduced by the Klein-Nishina effect, which will cause 
the Compton peak to shift to higher energies more slowly than 
the synchrotron peak (which depends on the square of the elec- 
tron energy). We then expect the TeV spectral shape to change 
quite little at the highest redshifts (z>0.1) - or equivalently, 
through the L>o.2Tcv(z) relation (see Fig.2e), at the highest 
luminosities (L>o.2TcV>10 45 erg cm" 2 s" 1 ). Conversely, for 
less luminous sources (i.e., sources caught in less active states) 
the energetics suggests a still largely Compton regime of the 
upscattered photons, which allows the Compton peak to shift 
to higher energies at the same rate as the synchrotron peak: 
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hence a spectral flattening is expected with increasing lumi- 
nosity for the more nearby sources. The data seem to conform 
to this prediction (see Fig.2f), with the weaker sources showing 
a steeper intrinsic value (o^ orr ~3) and the more powerful ones 
being spectrally revealed in the vicinity of the Compton peak 

The radio emission, on the other hand, is largely decoupled 
from the X-ray and TeV emissions. In the 'canonical' model of 
blazar emission the high-energy (X-ray and 7-ray) emission, 
though constituting most of the output, is produced in a region 
of the blazar's jet that has to be compact in order to account 
for the short-time variability shown by the blazars. This makes 
the radio emission from such region strongly self-absorbed, 
so that the observed radio emission must originate in a much 
larger volume. Therefore the link between the radio emission 
and the X-Ay-ray region may be much more subtle and less di- 
rect than suggested by a simple application of the SSC model 
(e.g., Costamante & Ghisellini 2002). As a consequence, the 
above arguments linking TeV slope and X-Ay-ray luminosity 
do not apply in the radio case: indeed, there is no correlation 
between intrinsic TeV slopes and radio luminosities (Fig.2g). 
There is, however, an overall trend of steeper observed slopes 
with increasing luminosities (Fig.2h), a consequence of more 
luminous objects being typically more distant and hence more 
IBL-affected at TeV energies, coupled with the essentially flat 
distribution of a^ orr (L5GHz)- 

5. Discussion and Conclusion 

Within the quite modest statistics and redshift range spanned 
by known TeV blazars, we find - contrary to some early expec- 
tations - no obvious dependence of the observed TeV slopes on 
redshift. However we find that, once the TeV slopes and lumi- 
nosities are corrected for intergalactic absorption, flatter slopes 
tend to be found in more distant, more luminous blazars. (A 
correlation of flatness and luminosity with redshift has been 
suggested also by Stecker et al. 2007.) We argue that this latter 
correlation descends from a selection effect: use of X-ray data, 
both simultaneous with TeV data and archival (when available), 
suggest that, for 0.1, more distant TeV blazars are in more 
active states (at the epoch of their observation); and more ac- 
tive states mean - given the link between blazar activity and 
spectral shape observed in nearby objects (e.g., Mrk 501) - 
harder TeV spectra. This activity-distance link, in turn, results 
from the selective effect of IBL absorption (that affects TeV, 
not keV, photons), which forces more distant objects to have 
higher intrinsic TeV/X-ray ratios in order to be detectable. In 
sum, the IBL appears to affect the blazars' TeV slope versus 
redshift distribution in two ways: first, by causing only objects 
in progressively more active states, hence with flatter intrinsic 
slopes, to be detected at higher redshifts; and then, by mak- 
ing the observed slopes steeper to the observer, more so for 
higher redshifts. The result of this twofold action of the IBL is 
- within the relatively small redshift range probed, 0.25 - 
a substantial lack of correlation between measured TeV slopes 
and redshift. 

It is clear, however, that if higher-luminosity/activity 
sources are detected at higher redshifts, the underlying Klein- 



Nishima regime will imply a substantially flat slope, a c °"^2, 
over a large luminosity range - hence over a relatively large 
redshift interval. Due to the redshift dependence of the IBL 
opacity, this entails a sharp correlation between observed spec- 
tral slope and redshift. We then predict that IACT observations 
of z^, 0.2 blazars will observe progressively higher a 7 's with 
increasing redshift. 

Of course, the above considerations are only as good as 
the adopted IBL model (Stecker et al. 2006) is. For a dif- 
ferent IBL redshift dependence, the resulting distribution of 
points in Figs.2a,b would also be different. More exotic in- 
terpretations of the intergalactic extinction, e.g., in the frame- 
work of the mixing between photons and a very light axion-like 
particle (De Angelis et al. 2007), would predict more mod- 
erate flux dimming and spectral distortion than entailed by 
most IBL schemes (e.g., Kneiske 2004; Mazin & Raue 2007). 
Reversing the arguments, if the maximum allowed hardness of 
TeV blazar spectra is known, the observed spectrum of a source 
with known distance in principle allows one to deduce an upper 
limit to the IBL. Aharonian et al. (2006b), using the H.E.S.S. 
observation of 1ES 1 101-032 and arguing that a^ orr >1.5, place 
an upper limit to the IBL at 1 .5^im that is close to its lower limit 
(i.e., the integration of galaxy counts), hence suggesting a 'low' 
IBL. Stecker et al. (2007) however, argue that o^ orr < 1 .5 can be 
produced if particles in blazar jets are accelerated at relativis- 
tic shocks, and hence loosen Aharonian et al.'s (2006b) upper 
limit, so permitting a higher IBL density field. Finally, if both 
the intrinsic blazar spectrum and the IBL are assumed, the ob- 
served TeV spectrum can be used to estimate the distance to 
the source (e.g., Mazin & Goebel 2007). 

It will be important to increase the statics in order to clarify 
the situation. This can be achieved both by probing larger red- 
shifts with deeper Cherenkov observations and by measuring 
TeV fluxes and spectra of nearby (z<0.1) blazars in different 
states of activity. The latter point, in particular, will enhance 
the statistics of the low-z portion of Figs.l,2a,2b,2f by adding 
crucial information on the quiet states of (nearby) blazars. As 
for the former point, the newly operating VERITAS telescope 
and the upcoming enhanced MAGIC and H.E.S.S. telescopes 
should improve the sensitivity of ground-based observations 
considerably, so adding to the number of known z>0.2 TeV 
blazars with spectral information above ^0.1 TeV. The im- 
pending GLAST satellite, too, will have an importat role in 
such expansion: its (largely IBL-unaffected) operational band, 
~0.01-100 GeV (with a sensitivity higher, by a factor of ^50, 
than its predecessor EGRET), will sample the blazar SED 
around the Compton peak, so in principle allowing the solu- 
tion of the SSC model and hence the knowledge of the intrinsic 
slope above ^0.1 TeV. 

Finally, we emphasize the need for future TeV blazar dis- 
coveries to be immediately followed-up in X-rays, that will 
help to validate or invalidate the main argument made in this 
paper, i.e. that more flaring sources are detected at higher red- 
shifts. This, if confirmed, would lend further support for the 
SSC process in blazar emission. 
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